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I.  INTRODUCTION 


In  recent  years  an  effort  has  been  made  to  determine  the  collision  cross 
sections  for  the  vibrational  excitation  of  1^0  by  0(^P).  This  work  has 
consisted  of  both  experimental  measurements  and  theoretical  calculations.  The 
measurements,  which  were  carried  out  by  Dunn  and  co-workers  in  a  rather  diffi¬ 
cult  shock  tube  experiment,  could  provide  only  a  limited  amount  of  data.^  In 
order  to  extend  these  data,  a  theoretical  program  was  initiated  by  Redmon  and 

r>  / 

Bartlett.  ~  They  calculated  an  accurate  O-f^O  potential  energy  surface 
and  then  used  this  surface  to  calculate  the  inelastic  cross  sections  by  the 
quasiclassical  trajectory  method.  The  results  they  obtained  were  in  fair 
agreement  with  the  limited  experimental  data  available.  This  is  encouraging; 
however,  it  is  not  definitive  since  both  the  theory  and  experiment  are  subject 
to  a  great  deal  of  uncertainty.  Other  theoretical  calculations  and  experi¬ 
mental  measurements  would  be  very  useful  in  judging  their  ultimate  accuracy. 
For  this  reason  we  have  decided  that  an  approximate  quantum  mechanical 
calculation  of  the  0  +  ^0  collision  problem,  employing  the  same  potential 
surface  as  the  classical  trajectory  calculation,  would  be  a  useful  test.  In 
addition,  from  a  purely  theoretical  point  of  view,  we  are  comparing  two  widely 
used  techniques  and  testing  them  (against  each  other)  in  new  physical  regimes. 

In  this  report,  we  present  the  results  of  a  quantum  mechanical  calcu¬ 
lation  of  the  vibrational  excitation  cross  sections  for  the  reaction 

0(3P)  +  H20  (000)  -*•  H20  (n  n  n3)  +  0(3P) 

where  nj ,  n2 ,  and  n^  are  the  symmetric  stretch,  bend,  and  asymmetric  stretch 
vibrational  quantum  numbers,  respectively.  The  relative  collision  energies 
considered  in  this  study  range  from  0.5  to  3.0  eV,  and  the  final  vibrational 
states  are  (010),  (100),  and  (001).  We  use  the  vibrational  close-coupling- 
infinite  order  sudden  (VCC-IOS)  method  to  treat  the  collision  dynamics  and 
use  the  same  potential  surface  as  Redmon  and  Bartlett  in  their  classical 
trajectory  calculation.  Thus,  our  study  will  give  a  direct  comparison  of 


the  VCC-IOS  and  classical  trajectory  methods  for  this  important  problem. 

The  application  of  the  VCC-IOS  method  to  atomic  scattering  from  both 
linear^  and  nonlinear^  triatoraic  targets  has  been  pioneered  in  recent  years  by 
D.C.  Clary.  The  method,  which  was  originally  developed  to  study  atom  diatomic 
molecule  scattering,  incorporates  an  exact  close-coupling  treatment  of  the 
vibrational  degrees  of  freedom  along  with  an  IOS  approximation  treatment  of 
the  rotational  degrees  of  freedom.  This  method  is  justified,  theoretically, 
when  the  rotational  energy  spacings  are  very  small  in  relation  to  the  trans¬ 
lational  and  vibrational  energies,  and  has  been  applied  with  great  success  in 
such  cases.  When  the  molecule  is  a  hydride,  as  in  the  present  case,  the 
rotational  spacings  are  much  larger,  so  the  method  has  less  justification. 

This  pessimistic  evaluation  would  certainly  apply  if  we  were  attempting  to 
calculate  transitions  to  a  specific  vibration-rotation  state.  However,  in  the 
present  problem,  we  are  calculating  vibrational  transitions  without  regard  to 
the  final  rotation  state,  i.e.,  we  sum  over  the  final  rotations.  Empirical 
evidence  cited  subsequently  indicates  that  these  summed  cross  sections  can 
indeed  be  calculated  accurately,  even  for  hydride  molecules. 

Our  justification  for  applying  the  VCC-IOS  method  to  the  0-H20  problem  is 
based  on  the  success  this  method  has  had  in  describing  atom-diatomic  molecule 
and  also  atom-symmetric  top  inelastic  scattering  problems.  In  particular  we 
refer  to  the  two  rather  severe  tests  of  the  VCC-IOS  method  given  in  references 
7  and  8,  where  the  collision  energies  were  low  and  the  rotational  energy  level 
spacings  of  the  target  molecules — H2  in  reference  7  and  NH3  and  H2CO  in 
reference  8 — are  rather  widely  spaced  as  they  are  for  the  H20  molecule.  In 
reference  8  it  was  found  that  all  the  IOS  cross  sections  to  individual 
rotational  states  were  in  qualitative  accord  with  accurate  values,  and  for 
most  the  agreement  was  within  20%.  The  most  important  result  obtained  in 
reference  7  is  the  observation  that  the  IOS  cross  sections,  summed  over  final 
states,  are  accurate,  even  for  very  low  collision  energies.  If  we  sum  the 
cross  sections  given  in  reference  8  over  final  states,  we  can  verify  that  in 
this  case  also,  the  summed  IOS  cross  sections  are  quite  accurate. 


Another  justification  for  using  the  VCC-IOS  approximation  in  the  present 
problem  is  based  on  the  good  agreement  between  exact  classical  trajectory  and 
classical  sudden  approximation  calculations  of  vibrationally  inelastic  cross 
sections  for  the  O-I^O  collision.  Redmon  et  al.  report  that  the  classical 
exact  and  classical  sudden  approximation  results  agree  within  a  factor  of  2 
for  this  problem. 


II.  CALCULATIONS 


A.  POTENTIAL  ENERGY  FUNCTION 

The  O-H2O  potential  energy  surface  (PES)  used  in  this  investigation  is 
the  same  PES  employed  by  Redmon  and  Bartlett  et  al.^>3  to  carry  out  their 
classical  trajectory  study.  This  surface  is  an  analytic  fit  to  ab-initio 
generated  points,  which  they  calculated  by  means  of  many  body  perturbation 
theory.  They  made  several  fits  to  the  PES.  The  particular  fit  used  in  this 
study  is  the  surface  that  they  have  designated  as  E^a.  A  complete  discussion 
of  this  surface,  the  ab-initio  calculations,  and  the  fitting  procedure  is 
contained  in  the  reports  by  Redmon  et  al.^~^  and  references  therein. 

B.  DYNAMICS 

Since  the  application  of  the  VCC-IOS  method  to  atom-nonlinear  (X-ABC) 
triatomic  molecule  collisions  is  described  in  detail  elsewhere, ^  we  will 
mention  here  only  a  few  salient  facts  in  order  to  fix  the  notation  of  the 
present  report.  An  appropriate  set  of  coordinates  to  describe  this  problem 
are  the  three  Eulerian  angles  a,  0,  and  y,  which  determine  the  orientation  of 
the  Eckart  coordinate  frame  for  the  triatomic  molecule;  the  three  normal  mode 
vibration  coordinates  QL ,  Q2,  and  Q3 ,  which  determine  the  position  vectors  of 
the  three  molecular  atoms  (A,B,C)  in  the  Eckart  frame;^  and  finally  the  three 
spherical  coordinates  R,  9,  and  $  of  the  colliding  atom  (X)  relative  to  the 
Eckart  frame.  The  molecule  lies  in  the  xy  plane  with  its  center  of  mass  at 
the  origin  of  the  Eckart  frame.  The  Schroedinger  equation  in  this  coordinate 
system  is 


(- 


2  viR 


R  + 


3R 


2pR 


+  HABC  +  Vlnf'E)  ♦  =  0 


int 


(1) 


where  p  is  the  reduced  mass,  L  is  the  angular  momentum  operator  for  the 
incident  atom,  HAg£  is  the  Hamiltonian  for  the  isolated  triatomic  molecule,  E 
is  the  total  energy,  and  Vint  (which  is  a  function  of  Qj ,  Q2 ,  Q3 ,  R,  0,  and  $) 
is  the  interaction  potential. 


In  the  VCC-IOS  approximation,  the  vibration-rotation  Hamiltonian  H^g^  is 

10  v 

replaced  by  the  simpler  zero  angular  momentum  vibrational  Hamiltonian  ^gC’ 
and  the  angular  momentum  operator  is  replaced  by  the  eigenvalue  form 

+  1)  .  The  result  is  an  extreme  simplification.  All  the  angular 
variables  ot,  0,  y,  0,  and  $  are  decoupled  from  the  dynamics, 
although  9  and  $  remain  in  the  wavefunction  as  parameters.  The  resultant 
Schroedinger  equation  is  then  solved  by  the  usual  close-coupling  technique  in 
which  the  wavefunction  is  expanded  in  the  following  basis  function 
representation 

*  -  2  f J  (R|M)  rn  CQ1  »Q2 ,  Q3)/R  (2) 

n 

The  basis  functions,  (Qj,  Q2*  Q3).  are  the  zero  angular  momentum  vibration 
states  of  the  isolated  triatoraic  molecule,  i.e.,  they  are  the  eigenfunction 
solutions  of  the  equation 

Ksc  -  rn(V  V  «3)  ■  °-  <3) 

The  index  n  in  the  above  expressions  represents  the  collection  of  three 
normal  mode  vibrational  quantum  numbers  (n^n2n3).  Equation  (3)  is  solved  by 
the  method  described  by  Whithead  and  Handy^  in  which  the  basis  functions 
^n  ^1’  ^2*  ^3^  are  exPaiu*ed  in  a  "primitive  basis  set” 


1  2 
-  I  I 

Ji  *  0  j2  -  0 


N, 

3  j , 

l  a(njj  ,j2,j  )  Q 

J3  -  0 


*  exp  +  w2Q2  +  u>3Q^)/2j  (4) 

Here  the  are  integers,  and  are  the  fundamental  normal  mode  vibration 
frequencies  of  the  triatomic  molecule. 
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Equation  2  is  substituted  into  the  VCC-IOS  simplified  version  of  Eq.  (1), 
multiplied  by  ?n  (Q^,  ,  Qj).  and  integrated  over  the  {Q^J  to  obtain  the 

vibrational  close-coupled  differential  equations 


where 


and 


[1^  +  k2  _  AU±il]  f*  (R;e,*) 

dRZ  R 


N 


H r  E  v  CR; 0 , 4>)  ff,(R;0,<fr) 

*2  n'=l  n,n  n 


(5) 


k2  -  ^  ( E-E  ) 

4t2 


(6) 


/  /  /  dQ1dQ2dQ3  Q2,  Q3) 


WVW  R’  9’  ^  W  V  V  (7) 

The  decoupled  variables  £,  0,  and  <f»  appear  as  parameters  in  Eq.  (5).  For 

each  assigned  set  of  parameter  values,  the  coupled  equations  are  integrated, 

subject  to  standard  scattering  boundary  conditions,  to  calculate  the  S  matrix 
g 

elements,  Sn  ^,(0,  $)•  Finally,  the  vibrational  excitation  cross  section, 
averaged  over  initial  rotation  states  and  summed  over  the  final  rotation 
states,  is  given  by  the  rather  simple  expression^ 


1  it/2  i  7 

[  f  dxd<p  Z  ( 2£+l )  I  <5  ,  -  S  ,(x,<f>)r 


a(n+n') 


(8) 


where  x  *  cos  0.  We  have  made  use  of  the  symmetry  of  1^0  molecule  to  reduce 
the  integration  range  of  both  the  x  and  4)  variables. 

C.  NUMERICAL  TECHNIQUES 

The  computer  code  used  to  carry  out  the  calculations  is  basically  the 
same  as  that  used  by  Clary**  in  his  He  +  SC>2  calculations  except  for  the 
modifications  described  subsequently. 

We  determined  that  it  was  not  necessary  to  integrate  the  scattering 

equations  for  every  £  value  in  the  range  0  <  £  <  £  .  Instead,  the  code  was 

max 

modified  to  calculate  for  every  tenth  £  value  and  to  interpolate  the  inter¬ 
vening  values.  This  resulted  in  a  substantial  savings  in  computation  time. 

The  numerical  quadrature  in  Eq.  (8)  was  changed  from  a  two-dimensional 
Simpson  rule  to  a  Monte  Carlo  procedure.  Using  this  method,  one  can  easily 
add  additional  integration  points,  if  needed,  without  discarding  previously 
calculated  points.  This  consideration,  combined  with  our  modest  accuracy 
requirement  that  the  integral  be  converged  to  about  20%  made  the  low  order 
Monte  Carlo  procedure  a  practical  alternative  to  the  higher  order  Simpson 
rule.  Using  this  method,  Eq.  (8)  becomes 


,  M  I* 

«(“■*')-  “4  I  I  (2*+l)l6nn’-Snn’  (xi'*i)l  (9 

Mlt  i-1  £ 

n 

where  the  M  values  of  x^  and  are  selected  at  random  with  uniform  dis¬ 
tribution  in  the  ranges  0  <  xt  <  1  and  -  tt/2  <  ^  <  ir/2.  We  selected  M  »  25, 
which  gives  a  nominal  convergence  of  20%. 

D.  CONVERGENCE  TESTS 

A  number  of  input  parameters  to  the  computer  code,  which  control  the 
accuracy  of  the  calculation,  must  be  determined,  primarily  by  trial  and 
error.  As  a  rule,  any  change  in  parameter  values  that  increases  the  accuracy 
of  the  computed  cross  sections  also  increases  the  time  and  therefore  the  cost 
of  the  calculation.  In  general,  we  adjusted  the  following  parameters  to 
obtain  about  5%  convergence:  The  lower  and  upper  integration  limits  on  the 


radial  variable  R  in  Eq.  (5)  were  set  at  2  and  11  bohr,  respectively.  The 

actual  method  for  obtaining  the  scattering  solution  to  this  equation  was  by 

R-matrix  propagation.  To  carry  out  this  calculation,  the  integration  range 

was  divided  into  100  equally  spaced  propagation  sectors.  The  maximum  angular 

momentum  quantum  number  was  set  to  £  -  200.  This  value  was  determined  to 

max 

be  sufficient  at  Che  highest  relative  translational  kinetic  energy,  3.0  eV, 
and  was  used  in  most  of  the  calculations. 

A  number  of  trial  calculations  were  performed  to  test  for  convergence 
with  respect  to  the  basis  set  expansion.  In  these  tests,  we  varied  the  number 
of  primitive  basis  functions  by  varying  the  parameters  Nj ,  N2 ,  and  N3  [see 
Eq.  (4)]  and  also  by  varying  N,  the  number  of  channels  in  the  close-coupled 
equations  [see  Eq.  (5)].  We  have  determined  that  the  minimum  basis  set 
required  to  obtain  reasonably  well  converged  cross  sections  is  given  by 
Ni  *  N3  »  2,  N2  ■  3,  and  N  ■  10.  It  is  of  some  interest  to  compare  the  energy 
levels  of  the  10  basis  functions,  calculated  using  this  primitive  basis  set, 
to  the  accurate  values  of  the  vibrational  energy  levels  for  H20.  This  com¬ 
parison  is  given  In  Table  1,  where  the  states  are  labeled  by  the  standard 
vibrational  quantum  number  designation.  Although  the  calculated  levels  are 
not  accurate  by  spectroscopic  standards,  they  are  "good  enough,”  as  shown  by 
our  convergence  tests,  to  serve  as  a  basis  set  for  the  scattering  calcula¬ 
tions.  The  vibrational  quantum  numbers  nj,  n2,  and  n^  have  a  simple 
interpretation^ — they  count  the  number  of  nodes  in  the  wavefunction;  n^ 
counts  nodes  along  the  Qj  coordinate,  n2  counts  nodes  along  Q2,  and  03  counts 
nodes  along  Q3.  The  primitive  basis  set  we  are  using  has  a  maximum  of  three 
bend  nodes.  This  is  the  reason  the  0,4,0  state  is  missing  from  our  computed 
basis  functions. 


In  order  to  adequately  test  the  convergence  of  the  scattering 
calculations,  it  is  not  sufficient  to  simply  increase  the  number  of  channels, 
N,  but  we  must  also  ensure  that  in  doing  so  we  have  increased  the  number  of 
nodes  in  the  basis  set  for  each  of  the  three  vibrational  modes  that  is  being 
tested.  In  Figs.  1  and  2  we  show  the  changes  that  occur  in  the  cross  sections 
for  the  excitation  of  the  010  and  001  (i.e.,  the  bend  and  asymmetric  stretch) 
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Table  1.  H2O  Vibrational  Energy  Levels  in  eV 


(nl 

n* 

t 

,  n3) 

Calculated3 

Accurate'3 

0 

0 

0 

0 

. 

0 

0.  0 

0 

1 

0 

0.202 

0.198 

0 

2 

0 

0.401 

0.391 

1 

0 

0 

0.464 

0.453 

0 

0 

1 

0.478 

0.466 

0 

3 

0 

0.614 

0.579 

1 

1 

0 

0.664 

0.648 

0 

1 

1 

0.677 

0.661 

0 

4 

0 

- 

0.762 

1 

2 

0 

0.857 

0.838 

0 

2 

1 

0.874 

0.851 

aThese  levels  are  the  eigenvalues  obtained  by  solving  Eq.  (3)  using  the 
primitive  basis,  Eq.  (4),  with  Nj  ■  2,  N3  ■  2,  and  N2  *  3.  The  corresponding 
eigenfunctions  form  the  basis  set  used  in  the  close-couping  scattering 
calculations. 

^Energy  levels  of  sufficient  accuracy  were  calculated  using  Eq.  (11,268)  and 
Eq.(lll,57)  in  reference  13. 
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Fig.  1.  Plot  Showing  the  Convergence  of  the  010  Cross  Sections  with  Respect  to 
the  Expansion  Basis  Set.  The  Labels  are  explained  in  the  text. 


states  when  the  basis  set  is  enlarged  to  include  an  extra  node  in  the  respec¬ 
tive  vibrational  states.  In  these  figures,  the  cross  section  curves  are 
labeled  by  (Np^.NjjN),  which  completely  describes  the  basis  set  used  in 
its  calculation.  We  have  also  indicated  the  maximum  number  of  nodes  in  the 
expansion  basis  sets.  In  Fig.  1  this  refers  to  the  bend  nodes  along  C^,  and 
in  Fig.  2  it  refers  to  the  asymmetric  stretch  nodes  along  Q3.  The  cross 
sections  shown  in  Figs.  1  and  2  were  calculated  for  a  single  value  of  the 
relative  orientation  angles  9  and  (The  values  used  were  0  *  90  deg  and 
$  *  -  7  deg.)  We  have  not  shown  the  convergence  test  of  the  symmetric  stretch 
100;  it  is  very  similar  to  the  asymmetric  001  case. 

In  addition  to  showing  the  change  due  to  Increasing  the  number  of  nodes 
along  Q2  in  Fig.  1,  we  also  show  Che  large  change  in  cross  section  for 
excitation  of  the  bend  state  when  we  reduce  the  number  of  bend  nodes  from 
three  to  two.  This  clearly  shows  that  at  least  three  bend  nodes  are  required 
for  reasonable  convergence  of  the  010  bend  excitation  cross  section. 

Many  additional  tests  were  also  carried  out  in  which  we  altered  the 
primitive  basis  so  that  we  could  add  even  more  symmetric  and  asymmetric 
stretch  nodes  (at  the  expense  of  fewer  bend  nodes).  The  results  obtained 
indicate  that  the  symmetric  and  asymmetric  stretch  excitation  cross  sections 
to  the  states  100  and  001  are  essentially  insensitive  to  the  addition  of  these 
extra  nodes  to  the  basis.  This  test  is  somewhat  artificial,  since  we  had  to 
reduce  the  bend  states  in  the  basis  set  to  less  than  optimum;  however,  we 
believe  it  is  still  a  valid  test  of  convergence  of  the  stretch  excitation 
cross  sections  with  respect  to  the  addition  of  extra  stretch  nodes  to  the 
basis  set. 

Our  conclusion  is  that  we  believe  our  cross  section  calculations  are 
adequately  converged  with  respect  to  the  expansion  basis  set  and  that  the 
degree  of  this  convergence  can  be  judged  approximately  by  a  comparison  of  the 
curves  plotted  in  Figs.  1  and  2.  Using  these  criteria,  the  convergence  for 
the  bend  excitation  varies  approximately  from  about  2  to  20%,  and  the 
convergence  of  the  symmetric  and  asymmetric  stretch  varies  from  about  10  to 
80%,  as  the  relative  translational  energy  varies  from  1.0  to  3.0  <iV. 


Tests  were  also  carried  out  to  verify  convergence  of  the  Monte  Carlo 
integration  in  Eq.  (9).  The  number  of  integration  points  were  increased  to 
M  »  50  and  M  ■  100  for  each  of  the  two  collision  energies  E  *  2  and  E  »  3 
eV.  These  tests  verify  our  estimate  that  the  calculations  carried  out  using 
M  ■  25  have  converged  to  a  nominal  20%  relative  accuracy. 


III.  RESULTS  AND  DISCUSSION 


The  cross  sections  for  transitions  from  the  ground  state  to  each  of  the 

three  vibrational  states  010,  100,  and  001  of  1^0  are  shown  in  Fig.  3.  Also 
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shown  for  comparison  are  Redmon's  classical  trajectory  results  *  and  two 
experimentally  determined  points.^  Based  on  our  discussion  of  the  convergence 
with  respect  to  basis  set  expansion  and  also  the  convergence  of  the  Monte 
Carlo  integration  over  the  angles  9  and  $,  we  estimate  a  relative  error  in  our 
calculated  results,  which  varies  from  about  20  to  40%  for  the  010  curve  and 
from  30  to  100%  for  the  100  and  001  curves. 

The  low  values  of  these  error  estimates  apply  at  the  low  collision 
energies  and  increase  to  upper  values  at  3.0  eV.  To  these  errors,  we  must 
also  add  the  error  resulting  from  the  use  of  the  rotational  sudden 
approximation.  This  error  is  much  more  difficult  to  estimate.  If  we  assume 
that  the  sudden  approximation  is  correct  to  within  a  factor  of  2,  as  has  been 
reported  for  the  classical  calculation  of  this  system,  then  the  maximum 
combined  uncertainty  in  our  cross  sections  will  be  about  a  factor  of  4.  Even 
if  we  assumed  a  much  larger  error  factor,  it  is  obvious  from  Fig.  3  that  a 
very  serious  disagreement  exists  between  our  VCC-IOS  calculated  cross  sections 
and  the  classical  trajectory  results.  Both  calculations  used  the  same 
potential  energy  surface,  eliminating  this  as  a  source  of  any  difference. 

Previously,  a  similar  comparison  between  classical  trajectory  and  VCC-IOS 
calculated  cross  sections  was  carried  out  for  the  O-CO2  collision. ^  In 
that  case,  contrary  to  the  present  O-H2O  problem,  the  agreement  between  the 
two  methods  was  fairly  good,  i.e.,  within  a  factor  of  2.^ 
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Fig.  3.  Cross  Sections  for  the  Collisional  Excitation  of  Ground  State  H2O  by 
0(^P)  to  the  Single  Quanta  Vibrational  States  010,  100,  and  001. 

Lower  three  curves  are  present  results,  computed  by  VCC-I0S  method. 
Upper  three  curves,  taken  from  reference  2,  were  computed  by  classical 
trajectory  method.  Dotted  curve,  taken  from  reference  14,  is  a  later 
and  more  accurate  classical  trajectory  calculation  of  the  001  cross 
sections.  The  points  A  and  G  are  experimentally  determined  values  of 
the  001  and  010  cross  sections  from  Ref.  1. 


IV.  CONCLUSIONS 


Quantum  mechanical  calculations  of  the  O-HjO  collision  excitation  cross 
sections  using  the  VCC-IOS  method  have  been  performed  in  the  relative  energy 
range  from  0.5  to  3.0  eV.  The  results,  which  are  shown  in  Fig.  3,  are  much 
smaller  than  the  classical  trajectory  calculated  cross  sections  obtained  using 
the  same  potential  energy  surface.  Both  of  these  methods  represent  state-of- 
the-art  techniques  for  chemical  dynamic  calculations,  and  both  are  widely 
used.  The  present  problem  has  attempted  to  push  both  methods  to  new  limits  of 
applicability.  The  large  disagreement  shows  that  one  of  the  methods  (or 
possibly  both)  has  been  pushed  too  far.  Further  work  should  be  carried  out  on 
this  or  similar  model  systems  to  resolve  the  discrepency. 
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LABORATORY  OP  ERA. T I ONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  Investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  Is  the  technical  staff's  wide-ranging  expertise  and  Its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 


Aerophy8lcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of-f ield-of-view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  infrared  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spacebome 
computers,  fault-tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications;  micro- 
wave  semiconductor  devices,  mlcrowave/mllllmeter  wave  measurements,  diagnos¬ 
tics  and  radiometry,  mlcrowave/mllllmeter  wave  thermonic  devices;  atomic  time 
and  frequency  standards;  antennas,  RJ  systems,  electromagnetic  propagation 
phenomena,  apace  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 


